Directed random walks on hierarchic trees with continuous branching: a 

renormahzation group approach 
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We investigate the directed random walk on hierarchic trees. Two cases are investigated: ran- 
dom variables on deterministic trees with a continuous branching, and random variables on the 
trees constructed trough the random branching process. We derive renormahzation group (partial 
differential) equations for the branching models with binomial, Poisson and compound Poisson dis- 
tributions of random variables on the links of tree. These renormaliation group equations are new 
class of reaction-diffusion equations in 1-dimension. 
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I. INTRODUCTION 

Models on hierarchic trees are rather popular in statis- 
tical physics [ij. Hierarchic tree models with random 
quenched disorder are especially popular object 

of research, they are connected with Random Energy 
Model. Due to special geometry, these models could be 
solved trough recursive equations Considering ran- 
dom walks (directed polymers) on trees, the recursive 
equation for these models have been approximated with 
Kolmogorov- Petrovsky-Piscounov (KPP)equation. 

The investigation of statistical physics models with 
quenched disorder on hierarchic trees with continuous 
branching are the subject of extensive studies. The mod- 
els with continuous branching are especially popular in 
financial market literature as well as in turbulence. 

In % has been considered hierarchic trees with contin- 
uous branching, when there are random variable with 
normal distribution on the tree. These version of mod- 
els later has been investigated in [13: [HI , [12! > applying to 
the problems of disordered systems and string theory. To 
solve hierarchic models with continuous branching, in [^ 
has been used a renormahzation group equation, first ob- 
tained in [T3I . This equation is similar to KPP equation, 
see [l2| for discussion. In the current article we will de- 
rive renormahzation group equation for the general case 
of distribution for random variables on the branches of 
hierarchic tree. 

We are going to investigate the random directed walks 
on hierarchic trees described trough partition function: 



(1) 



where jji are defined on the endpoints of the tree and /3 
is the inverse temperature. 
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The model is defined by considering a tree with branch- 
ing number q and K levels of hierarchy, therefore there 
are = endpoints of the tree. At any level of hierar- 
chy there are q links from every node, and on every link 
a random variable ei is defined. The variable yi entering 
the definition of the partition function Eq.(l) is associ- 
ated with a unique path going from the origin O to the 
i-th tree endpoint and amounts to summing up variables 
ci along this path. 



(2) 



One can calculate mean free energy In Z using the iden- 
tity i 



< InZ >= I ( — 
/o P 
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r'(i) 



Inpd < e *^ > 



(3) 



Thus we need to calculate 

gip,(3) = Gk{x) =< exphp^e^^^'] > |„, 

i 



(4) 



for the average < exp[— p^ - e~^^'] > |e, over the config- 
urations on the K-level hierarchic tree. 

Having the expression for g{p, f3), we can calculate the 
moments of Z~": 



dpp"gip, P) 



(5) 



In Ref. [1] an equivalent problem has been considered, 
and recursive relations have been introduced . Following 
[5j, we define 



Go(x) 



exp[— e 



(6) 



and then other K functions Gi{x) 



exp[— e 



-Px 



for the model on the 
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trees with I level of hierarchy, < I < K. We identify 
5(p,/3) = Gx(-lnp//3). 

Consider the 1-level hierarchic tree. It could be frac- 
tured into q trees with hierarchy level I — 1 each. The yi 
at the endpoints of 1-level tree can be derived from the 
yi of the {I — 1) level tree adding a random variable e. A 
simple consideration gives a recursive relation 0, 0| : 



Gi{x) ^[ de Gi-i{x + e)p{e)] 



(7) 



Considering recursive relations for 1 < I < K with the 
boundary condition by Eq. (6), we can calculate Gk- 

Eq.(7), derived in |2j]-|5|] is an exact recursive relation. 
Similar recursive equations have been considered in the 
research of real space renormalization approach to quan- 
tum disorder in d-dimensional space, see Eq.(32) in [20| . 
The latter is some approximation, while Eq.(7) is an ex- 
act relation. The more serious difference- our random 
variables e have independent distributions for different 
hierarchy levels I, while at quantum disorder point [20] 
there is some correlation of noise at different hierarchy 
levels. To get some analytical estimates, in [13] has been 
considered the limit of large branching number q (K in 
notation of similar to those used in @. used 

general distribution of random variables. 

The KPP like equation has been applied for investiga- 
tion of quenched disorder model. In the next section we 
will consider the g — ?> 1 limit of Eq.(7), deriving renor- 
malization group (KPP like) equation for the general dis- 
tribution of e. 



II. RANDOM WALKS ON A HIERARCHIC 
TREE 

A. Renormalization-group equation for Hierarchic 
trees with continuous branching 

While Eq.(7) is derived for integer q, we can consider 
the equation for any positive value of q as well. Let 
us consider the model, where the branching number is 
close to 1. We can now derive exact differential equation 
instead of iteration equation Eq. ([7]). We have K = 1/ Av 
levels of branching on our tree with q = e^^"" , where 
LAw < 1. We identify 

Gi{x) = G{x,v), 
I 



g{p^P) = G{-^) 



(8) 



To define a distribution for random variable e;, let us 
start with some random distribution p{x). In the next 
step we calculate the distribution for the sum of L ran- 
dom variables with such distributions. For our purposes 
it is convenient to use the following representation for p: 



1 /"^ 

p{x) ~ — I dh exp[0(ift,) — ihx] 



(9) 



Then for the e, a sum of L random variables x, we get 
the compose distribution e) just multiplying (p in the 
exponent Of Eq. (9) by L: 



1 

2^ 



dh exp[L0(i/i) — ihe] (10) 



To ensure the probability balance condition, there is a 
constraint 



0(0) - 

We can consider the Taylor expansion 
0(zA:) = ^fe/(ifc)' 



(11) 



(12) 
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According to our notations p{x) = p{\,x). A concrete 
case of the latter distribution with In cosh(h) has been 
used while considering the diluted REM '19]. Here we 
will consider another distributions, popular in turbulence 
and financial theory. 

We take p{L, yi) as a distribution of yi. The L in the 
exponent of Eq.(lO) gives correct scaling for the < e^^* >, 



<e^" >=exp(i0(/3)). 



(13) 



where is the number of end-points of hierarchic tree. 

As yi is a sum of 1/Av random variables e on our tree, 
we define the following distribution for the distribution 
of random variables on the links p{e). 



p{e)=p{LAv,e) = 
1 f°° 

— / dh exp[(p{ih)LAv — ihe] 
27r 7_oo 



(14) 



For the LAv <C 1 we expand in the exponent the terms 
with bi,l > 1, and derive: 

/ deG{x + €,v)p{e) = — / dh exp[biLdv — ihe] 

x{l + LdvJ2bi{ihy)G{x + €,v) ^ I d€[S{e ^ biLAv) 
1 f°° 

+LAv^bi{ih^)— / d/i exp(-i/ie)]G(a; + e,iOl5) 



l>2 



We assume a smooth behavior of G(x, v) at the the in- 
finity X — >■ ±cx), and ^ {dx)^^ |x=±oo — for n>l. 
Integrating by parts we derive: 



deG{x + €,v)p{e) = G(a;, u) + iAz;0( — )G(x, w)(16) 



Thus we obtain: 

G{x, v) [(G -I- LAv (j}{dx)G{x, v)]^+^^'' « 
G{x,v) + Av L{(t){dx)G{x,v) + G{x,v)\nG{x,v)) (17) 
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dG{x,v) 
dv 



4>{dx)G{x,v) +G{x,v)\xvG{x,v) (18) 



where < u < 1 plays the role of time of reaction- 
diffusion equation, and we should solve the equation with 
the initial distribution 



G(x,0) =cxp[-e-''^ 



(19) 



Due to proper choose of the scaling in the exponent of 
Eq.(14), there is no L dependence in Eq.(18). 

Contrary to the KPP equation in the case of models 
of Ref. [2,] , Eq. (fT8| is an exact equation. 

Solving Eq.(18) for the initial distribution by Eq.(19), 
we can calculate p{Z) using the inverse Fourier transfor- 
mation. 



B. Different distributions 

Consider different distributions, popular in cascade 
processes. 

The normal distribution. Taking 



we obtain 



(j){k) = fcV2 



dG{x,v) d^G{x,v) 



dv 



dx"^ 



GlnG 



(20) 



(21) 



We derived Eq. (21) in [9!] solving the continuous branch- 
ing model. First time this equation has been derived in 
[13| for investigation of non-linear diffusion. 
The binomial distribution. We have 



i^2<5(e - (1 — ai) ai5{e — {I + a2)) 



Oil + Oi2 

We get for (j){k) an expression: 



^(fc) = fc(l - ai) +ln[- 



Ot2 



(22) 



-ln[l H e 



E 



(-1)" ,a2 fc(ai-Ha2)l 



(23) 



Eventually we get the following PDE: 



^^("'")-GlnG+(l-aO-^^^"'"^ 



ln[ 



0.2 



dv dx 

J2^-^[—Giv,x + n{ai+a2)] (24) 



n>l 



We are looking solutions where G{x, u) — ^ for x — oo, 
therefore our equation is well defined. 



Gamma distribution. Now we consider the follow- 
ing distribution for positive e 



r(7) 



-7£ 



(25) 



We have the following expression for (t>{k): 

0(fc) =7ln(l + fc/7) (26) 
While we can formulate formally a PDE in this case, 

^^^=7ln[l + -^] G(a;,i;)+GlnG, (27) 
ov 7 ox 

the equation can be better formulated in Fourier space. 

The Poisson distribution. It is popular in financial 
mathematics literature. We have integer values of e 



/5(e) - e- V/e! 



(28) 



Eq. dH) gives (j){k) = 7(-l + e*^). Then from Eq. ^ we 
derive: 



dG{x,v) 
dv 



7[-G + G(a; + l)] + GlnG (29) 



This renormalization group equation differs from Eq. (j2ip 
in that the second order derivative is replaced by finite 
difference. 

Compound Poisson distribution. The cascade 
processes with such distribution are also rather popular 
in literature. 

Now we have Poisson distribution for integers n given 
by Eq.(28) i.e. pi{n) = e^'^^^/nl and define e as a sum 
of n random variables xi with some random distribution 
p{x), defined trough the representation: 

1 

p{x) = — / dh exp[a{ih) — ihe] (30) 
27r J_oo 



{k)=\n p{y)dy exp[ky] (31) 



Thus 



A simple calculation gives the function (j){k) for the dis- 
tribution of e 

(l){k) = 7(e"W - 1) (32) 
Thus we get the following PDE for G{v,x): 

=^[e"(7fe)-l]G(a;)]+GlnG (33) 

ov 

or, using Eq.(31), 
dG{x, V 



dv 



7[ / dyp{y)Gix + y,v)- G(x, v)] + G In (^4) 
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C. The phase transition point 

General case of distribution. The free energy of 
directed polymer model Eq.(7) for any q is equivalent 
to the free energy of corresponding REM with the same 
number of energy configurations and distribution oi m. 
This issue is well discussed in REM literature 0, H, [E, 0] ■ 

We have M = endpoints of the tree, and any of yi 
has a distribution by p{L,y), see Eq. (10). Let us con- 
sider REM with M energy level with independent distri- 
bution of energy by Eq.(lO). 

At high temperature phase, assuming < InZ >= In < 
Z >, we have fsj. 



KlnZ >= L{l + (j){l3)) 



(35) 



The phase transition is at the point where the entropy of 
free energy by Eq.(35) disappears 3, 



X/?e)-/3c0'(/3c)=O 



(36) 



Binomial distribution. Now we have the following 
expression for the free energy in the high temperature 
phase 



\nZ a.2 
—— = in[ ( 

L ai + a2 



,/3(l-ai) 



ai 



p/3(l+a2)l 



ai + a2 



(37) 



Eq.(36) has no solution in this case, therefore system is 
always in the high temperature phase. 

Gamma distribution. Now we have the following 
expression for the free energy in the high temperature 
phase 



InZ ^ ,^ /3, 
— = 1 + 7 
L 7 



(38) 



Eq.(36) has no solution in this case, therefore system is 
always in the high temperature phase. 

The Poisson distribution. For this case we derive 
for the free energy of high temperature phase and tran- 
sition point: 



l + 7(e^-l); 



\nZ 

l + 7[e^=(l-/3e)-l] = 



(39) 



In the next sections we will use the travelling wave anal- 
ysis to investigate the paramagnetic phase solution of 
our equations and identify the transition point. For in- 
vestigation of the spin-glass like solutions one needs in 
advanced mathematics of jT3|. 

Compound Poisson distribution. Now we have for 
the high temperature phase the following expression for 
the free energy 



^ = l+7(e"(«-l) = 
l + 7(/ dxp{x)e^'' - I) 



(40) 



For the critical point we get an equation: 

1 + 7[ / dxp{x)ef^'''' - 1 - /3 / dxp{x)ef^''''x\ = (41) 



D. 



KPP versus new renormalization group 
equation 



Let us try to describe our hierarchic tree model using 
KPP equation 



dQ{x,v) _ d^Qix,v) 



(42) 



dv dx^ 

In f3| has been derived Eq.(42) as an approximation of the 
Eq.(7) for the special case of Eq.(20). Later Eq.(42) has 
been derived in [T^ as an approximate renormalization 
group equation for the models of random disorder with 
logarithmic correlations. While in [1, [1] has been already 
considered the general distribution, the KPP has been 
derived only for the normal distribution. In po| , while us- 
ing general distribution of random variables, the authors 
does not used KPP. The reason is clear: for the general 
distribution of random variables on the tree, the KPP 
equation gives wrong results even for transition point, 
contrary our new renormalization group equation which 
is exact. 

To compare we take the case by Eqs.(28) with 7 = 1. 
The KPP version, having the same variance of distribu- 
tion p(e), gives 



dG{x,v) dG{x,v) d^G{x,v) 



dv dx 
Using the mapping: 



dx^ 



G{v,x) = Q{v,x + t) 



(43) 



(44) 



we return to the KPP equation by Eq.(42). This equation 
is well investigated in literature. It describes a transition 
between spin-glass and paramagnetic phases at 

/3c = \/2 (45) 

and free energy at paramagnetic phase is 

<lnZ >=L{l + l3 + l3^) (46) 

The correct transition point of the model and free energy 
are: 



/3c = 1 
< InZ >= Le^ 



(47) 



Thus the KPP wrongly describes the model, and we need 
in new exact renormalization group equation for the gen- 
eral form of distribution of random variables on the tree. 

Let us investigate our Eq.(29). G{v,x) is a monotonic 
function of x according to the definition by Eqs.(4),(7), 



G{v,x)^l, for a; — > 00, 
G(f,a;)— T^O, for x —00 



(48) 
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Consider the case of small p in Eqs.(4),(7). The direct 
calculations give 

G{v,x) = exp[-pexp/3x + + 7(e^ - 1))], (49) 

The latter expression is not a solution of Eq.(43). Nev- 
ertheless, it is the solution of Eq.(29) when 



pexpl^x + v{l + j{e'^ - 1))] < 1 



(50) 



Let us investigate the behavior of the asymptotic so- 
lution for the more general situation than the case of 

Eq-(50). 

Following to [2|,l5|,llJ], we assume a traveling wave like 
asymptotic solution 

G{v,x) = g{x + cv) (51) 

We get ODE: 

cg'ix) = [~gix) + g{x + 1)] + g{x) \ng{x) (52) 

Eq.(52) describes a traveling wave with a front at 

X = —cv (53) 

We assume that for large negative x there is a solution 

g{x) = exp[-e^T (54) 

Putting the latter expression to Eq.(52), we get the fol- 
lowing equation for c: 



cip) 



[e^ - 1]7 + 1 

/3 



(55) 



Eqs.(51),(54),(55) give the free energy expression by 
Eq.(46). 

We can identify the phase transition point as the value 
of (3 giving the maximum velocity, looking the equation: 



c'(/3c) = 



(56) 



Eq.(56) gives the result of Eq.(47), f3c ~ I for the case 
7 = 1. 

At high value of (3 there is another solution for c(/3), 
but for its derivation one needs in advanced mathemati- 
cal approach of |17| • 

Thus the investigation of our renormalization g roup 
equation, using the qualitative approach of 0, [3, Il4{ . 
gives exact expression for the free energy. 



III. RANDOM BRANCHING MODEL 

In the previous section we derived new renormalization 
group equation, using continuous deterministic branch- 
ing, later considering (7 — > 1 limit. Let us derive a sim- 
ilar renormalization group equation, using less abstract 
model of random branching 



We again have a tree. It starts at some point O, and 
grows down, where the vertical coordinate measures the 
time V. We put random variable e; with distribution by 
Eq.(8) on the links (between two adjacent nodes of the 
tree), replacing Av by the difference of time coordinate 
between two nodes of the branch. During period of time 
dv there is a branching with a probability dt. 

Consider the dynamics of the partition sum Z(/3,ti) 
defined as 



(57) 



Here the index j numerates the endpoints of our tree at 
the time v. 

Z{/3,v) has a deterministic behavior at v = 0: 



Z(/?,0) = 1 



(58) 



Z{/3, v) is a random variable at w > both due to the 
randomness of branching and randomness of e; . 

Then one has the following recursive equation for the 
random partition function , 

Z{(3,v + dv) = Z{/3,v)e-^\ (59) 

with probability 1 — dv, and 

Z{I3, v + dv)^ iZ\l3, v) + Z\/3, v))e-^' (60) 

with probability dv. Here Z^{P,v) and Z^(/3,u) are ran- 
dom independent variables with the same probability dis- 
tribution us Z(/3,u). 

We identify Z(0,f) with the number of endpoints of 
the tree, and Eqs.(59)-(60) give the intuitive result 



< Z{0,v) >= exp(i). 



(61) 



confirming the self-consistence of the choice Eqs.(59)- 
(61). Actually the random process Z{t) is defined com- 
pletely trough Eqs.(58)-(60), and we can work with these 
equations missing any reference to the branching trees. 
We define now, following 



G{v,x) =< exp[-e-''^Z(i;)] >, 
Then we have an equation 



(62) 



G{v + dv, a;) = (1 - dv) J dep{e)G{v, x + e) + 

dvG{v,xf (63) 
Eventually we get the following equation: 
dG{x,v) 



dv 



(j){dx)G{x, v) ~ G{x, v){l - G{x, v)) (64) 



For the binomial distribution we get: 

^^(^'")=-G(l-G) + (l-«0-'^(^'^) 



ln[ 



dv 

(\2 



dx 



ai + a2 



J2^^[—G{v,x + nia,+a2)] (65) 



ni ot\ 
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For the Poisson distribution we get 
dG{x,v) 



dv 



= 7hG + G{x + 1,)] + G{x, vY ~ G{x, w)(66) 



We can repeat the derivations of the section II-D, and 
again get the same solutions for the free energy of para- 
magnetic phase and the critical temperature. 
For the compound Poisson distribution we get: 



dG{x,v) 
dv 



l[ / dyp{y)G{x + y,v)- G{x, v)] 



+G(G-1) (67) 



IV. CONCLUSIONS 

We provided an exact (renormalization group) par- 
tial differential equation, for the analysis of models on 
hierarchic tree with continuous branching and general 
distribution of random variables on the tree. We de- 
rived new classes of reaction-diffusion equation Eq. (18), 
which sometimes has finite difference Eqs. (24), (29), (34). 
While recursive equations for the hierarchic models with 
general distributiqon of random variables on the tree have 
been well known 11],[2^, our article is the first result to 
replace theses recursive equation with correct PDF. 

Our derivation of renormalization-group equation is 
quite rigorous. The renormalization group equation is 
exact, including finite size corrections, while the renor- 
malization group equation in [l^| is an approximation 
and could be applied to calculate only bulk characteris- 
tics of corresponding models. 

The KPP equation, used in 0, [ij] , correctly (exact 
in the thermodynamic limit) describes the bulk free en- 
ergy and critical temperature for the hierarchic tree mod- 
els with normal distribution of random variable 0, Q, 



as well as for the finite dimensional models of disorder 
with logarithmic correlations [l3]- As we checked, KPP 
equation failed to describe correctly the hierarchic mod- 
els with non-normal distributions. Thats why the general 
case of hierarchic model, given by recursive equation (7) 
derived in 0, never has been investigated using KPP in 
d H. To verify our new renormalization group equa- 
tions, we used the idea of traveling wave and methods 
of [1, [i3]- Analyzing our equations, we found correct ex- 
pression for the free energy in paramagnetic phase and 
the phase transition temperature. 

An open problem is to construct the finite dimensional 
models of disorder, which could be described by our new 
equations, as the models of [111 are described by KPP. 
We hope to succeed using dynamic stochastic processes 
in 1 dimensional case. 

Our approach (continuous branching) could be applied 
to investigate the quantum disorder in d-dimensional 
space. In case of continuous branching we derived ex- 
act PDF equations which could be solved numerically, 
while in alternative approach of '2C?] has been derived 
approximate renormalization group equations and there 
are errors 0(1) in the expression of free energy (bulk term 

The normal distribution version Fq. (PT|) has a deep 
mathematical meaning: in some sense it describes the 
p-adic space with p = 1 [l^. We hope that the investi- 
gation of these new reaction-diffusion equations will be 
further encouraged. 

Our renormalization group equations could also be ap- 
plied to calculate also the complicated correlation func- 
tion, as has been done in Ref. [ll| for the string case. 
This work was supported by Academia Sinica, National 
Science Council in Taiwan with Grant Number NSC 99- 
2911-1-001-006, and National Center for Theoretical Sci- 
ences (Taipei Branch). 
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